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The magnon exchange mechanism of superconductivity was 
developed to explain in a natural way the fact that the su- 
perconductivity in UGe2, ZrZni and URhGe is confined to 
the ferromagnetic phase. The order parameter is a spin anti- 
parallel component of a spin-1 triplet with zero spin projec- 
tion. The transverse spin fluctuations are pair forming and 
the longitudinal ones are pair breaking. In the present pa- 
per, a superconducting solution, based on the magnon ex- 
change mechanism, is obtained which closely matches the 
experiments with ZrZn-i and URhGe. The onset of super- 
conductivity leads to the appearance of complicated Fermi 
surfaces in the spin up and spin down momentum distribu- 
tion functions. Each of them consist of two pieces, but they 
are simple-connected and can be made very small by vary- 
ing the microscopic parameters. As a result, it is obtained 
that the specific heat depends on the temperature linearly, at 
low temperature, and the coefficient 7 = ^ is smaller in the 
superconducting phase than in the ferromagnetic one. The 
absence of a quantum transition from ferromagnetism to fer- 
romagnetic superconductivity in a weak ferromagnets ZrZn^ 
and URhGe is explained accounting for the contribution of 
magnon self-interaction to the spin fluctuations' parameters. 
It is shown that in the presence of an external magnetic field 
the system undergoes a first order quantum phase transition. 
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Very recently ferromagnetic superconductivity (f- 
superconductivity) has been observed in UGe2 [Q, 
ZrZn2 Q and URhGe [||. The superconductivity is 
confined to the ferromagnetic phase. Ferromagnetism 
and superconductivity are believed to arise due to the 
same band electrons. The persistence of ferromagnetic 
order within the superconducting phase has been ascer- 
tained by neutron scattering. The specific heat anomaly 
associated with the superconducting transition in these 
materials appears to be absent. 

At ambient pressure UGe2 is an itinerant ferromag- 
net below the Curie temperature Tc — 52K, with low- 
temperature ordered moment of fig = lAfis/U. With 
increasing pressure the system passes through two suc- 
cessive quantum phase transition, from ferromagnetism 
to f-superconductivity at P ~ 10 kbar, and at higher 
pressure Pc ^ 16 kbar to paramagnetism ||l|,Q.At the 
pressure where the superconducting transition tempera- 
ture is a maximum Tgc — 0.8K, the ferromagnetic state is 
still stable with Tc = 32K, and an ordered moment about 
1.0^b/C/ The specific heat coefficient 7 = § 

increases steeply near 11 kbar and retains a large and 
nearly constant value 



The ferromagnets ZrZn^ and URhGe are supercon- 
ducting at ambient pressure with superconducting criti- 
cal temperatures Tsc = 0.29K and Tgc = 0.25K respec- 
tively. ZrZn2 is ferromagnetic below the Curie tem- 
perature Tc = 28.5if with low-temperature ordered mo- 
ment of /is — O.llfiB per formula unit, while for URhGe 
Tc = 9.5K and fig — 0.42/is. The low Curie temper- 
atures and small ordered moments indicate that com- 
pounds are close to a ferromagnetic quantum critical 
point. A large jump in the specific heat, at the temper- 
ature where the resistivity becomes zero, is observed in 
URhGe. At low temperature the specific heat coefficient 
7 is twice smaller than in the ferromagnetic phase. 

The most popular theory of f-superconductivity is 
based on the paramagnon exchange mechanism [^,0. 
The order parameters are spin parallel components of the 
spin triplet. The superconductivity in ZrZn2 was pre- 
dicted, but the theory meets many difficulties. In order 
to explain the absence of superconductivity in param- 
agnetic phase it was accounted for the magnon param- 
agnon interaction and proved that the critical temper- 
ature is much higher in the ferromagnetic phase than 
in the paramagnetic one To the same purpose, the 
Ginzburg-Landau mean-field theory was modified with 
an exchange-type interaction between the magnetic mo- 
ments of triplet-state Cooper pairs and the ferromagnetic 
magnetization density |^.In ||l^ the authors make the 
important assumption that only majority spin fermions 
form pairs. Then, only minority spin fermions contribute 
to the asymptotic of the specific heat, and the coefficient 
7 = ^ is twice smaller in the superconducting phase. The 
result closely matches the experiments with U RhGe [|j , 
but does not resemble the experimental results for U Ge2 
and ZrZn2- The assumption seems to be doubtful for 
systems with very small ordered moment. Despite of the 
efforts, the improved theory of paramagnon induced su- 
perconductivity can not cover the whole variety of prop- 
erties of f-superconductivity. 

In the present paper an itinerant system is considered 
in which the spin-^ fermions c„{x){c^{xy) responsible 
for the ferromagnetism are the same quasiparticles which 
form the Cooper pairs. The exchange of spin fluctuations 
leads to an effective four fermion theory. It describes 
the interaction of the components of spin-1 composite 
fields (tT7Ti + itiii) which have a projection of spin 
1,0 and -1 respectively, and the interaction of the spin 
singlet composite fields ti — it- The spin singlet fields' 
interaction is repulsive and does not contribute to the 
superconductivity |jl^. The spin parallel fields' interac- 
tions are due to the exchange of paramagnons and do 
not contribute to the magnon-mediated superconductiv- 



ity. The relevant interaction is that of the ti + it fields. 
The potential of this interaction has an attracting part 
due to exchange of magnons and a repulsive part due to 
exchange of paramagnons. 

By means of the Hubbard-Stratanovich transformation 
one introduces ti + it composite field and then the 
fermions can be integrated out. The obtained free en- 
ergy is a function of the composite field and the integral 
over the composite field can be performed approximately 
by means of the steepest descend method. To this end 
one sets the first derivative of the free energy with respect 
to composite field equal to zero, this is the gap equation, 
and looks for a solution which minimizes the free energy. 

The gap is an antisymmetric function A(— fc) = 
— A(fc), so that the expansion in terms of spherical har- 
monics Yim(Xlj:) contains only terms with odd I. I assume 
that the component with I = 1 and m = is nonzero and 
the other ones are zero 



A(fc) = Aio(fc)W — cos( 
An 



(1) 



Expending the potential in terms of Legendre polynomial 
Pi one obtains that only the component with I = 1 con- 
tributes the gap equation. The potential Vi {p, k) has the 
form, 
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where M is zero temperature dimensionless magnetiza- 
tion of the system per lattice site and p is the spin stiff- 
ness constant which is proportional to M [p — Mpo) The 
constants /3, po Siiid b are phcnomcnological ones subject 
to the relation (3 = — ^ > ^^^^ ~ i ^' 
where the parameter r is the inverse static longitudi- 
nal magnetic susceptibility, which measures the deviation 
from quantum critical point. A straightforward analy- 
sis shows that for a fixed p , the potential is positive 
when k runs an interval around p (p — A,p + A), where 
A is approximately independent on p. In order to al- 
low for an explicit analytic solution, I introduce further 
simplifying assumptions by neglecting the dependence of 
Aio(fc) on k (Aio(fc) = Aio(p/) = A,p/ = ^/2pm) and 
setting Vi{pf,k) equal to a constant Vi within interval 
{pf — A,pf + A) and to zero elsewhere. 

To ensure that the fermions which form Cooper pairs 
are the same as those responsible for spontaneous mag- 
netization, one has to consider the equation for the mag- 
netization 
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as well. Then the system of equations for the gap and for 
the magnetization determines the phase where the super- 
conductivity and the ferromagnetism coexist. The system 



can be written in terms of Bogoliubov excitations, which 
have the following dispersions relations: 
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where A(fc) is the gap J is the spin exchange con- 
stant, and e(fc) = ^ — P- 

At zero temperature the equations take the form 
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where t — cos 0. 

The solution of the system which satisfies ^^A < JM 

is discussed in ||l2[. In the present paper one looks for a 
solution of the system which satisfies 



-A > JM 
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One is primarily interested in determining at what 
magnetization a superconductivity exists. The inequal- 
ity Eq. (^) shows that the gap can not be arbitrarily small 
when the magnetization is finite. Hence the system un- 
dergoes the quantum phase transition from ferromag- 
netism to f-superconductivity with a jump. Approach- 
ing the quantum critical point from the ferromagnetic 
side, one sets the gap equal to zero in the equation for 
the magnetization and considers the gap equation (|^) 
with magnetization as a parameter. It is more convenient 
to consider the free energy as a function of the gap for the 
different values of the parameter M . To this purpose I 
introduce the dimensionless "gap" x and the parameters 
s, A and g 
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Then the free energy is a function of x and depends on 
the parameters s, A and g. 
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The dimensionless free energy F(x) is depicted in Fig.l 
for A ~ 0.08, g — 20 and three values of the parame- 
ter s, s — 0.8, s = 0.69 and Scr = 0.595. As the graph 



shows, for some values of the microscopic parameters A 
and g, and decreasing the parameter s (the magnetiza- 
tion), the system passes trough a first order quantum 
phase transition. The critical values Scr and Xcr satisfy 
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FIG. 1. The dimensionless free energy F{x) as a function 
of dimensionless gap x. \ = 0.08, g = 20, Si = 0.8(upper 
line), S2 = 0.69(middle line) and Scr ~ 0.595(lower line). 



Eqs (|ri|,^2|) are the solution of the system Eqs.(||j6|) near 
the quantum transition to paramagnetism. The second 
derivative of the free energy Eq.(||) with respect to the 

gap is positive when > (^[f-)^, hence the state 

where the superconductivity and the ferromagnetism co- 
exist is stable. 

When superconductivity and ferromagnetism coex- 
ist, the momentum distribution functions n^{p,t) and 

{p, t) of the spin-up and spin-down quasiparticles have 
complicated Fermi surfaces. One can write them in terms 
of the distribution functions of the Bogoliubov fermions 

(p, t) = (p, t)m (p, t) + (p, t)n2 {p,t) (13) 
uHp, t) = u^{p, t){l - ni{p, t j) + v^{p, t){\ - n2{p, t)) 

where it(p, t) and u(p, t) are the coefficients in the Bogoli- 
ubov transformation. At zero temperature ni{p,t) = 1, 
n2{p,t) = Q{—E2{p,t)), and the Fermi surface Eq.(|l0|) 
manifests itself both in the spin up and spin-down mo- 
mentum distribution functions. The functions are de- 
picted in Fig. 2 and Fig. 3. 



Varying the microscopic parameters beyond the critical 
values, one has to solve the system of equations (||||).The 
equation of magnetization (|5|) shows that it is convenient 
to represent the gap in the form A = y/^K{AI)JM , 
where k{M) > 1. Then the equation E2{k,t) = 0, which 
defines the Fermi surface, has no solution if — 1 < t < 



k(M) 



and 



k(M) 



< i < 1, and has two solutions 
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when 



< t < 



Kp?) ^ " ^ K{M y "^^"^ solutions (|T^) determine 
the two pieces of the Fermi surface. They stick together 
sX t = ±-j^pjj, so that the Fermi surface is simple con- 
nected. The domain between pieces contributes to the 
magnetization M in Eq.(|^), but it is cut out from the 
domain of integration in the gap equation Eq.(|6|). 

When the magnetization approaches zero, one can ap- 
proximate the equation for magnetization Eq.(^) substi- 
tuting from Eq.(|lO|) in the the difference (Py )^ — (pj)^ 
and setting = Pf elsewhere. Then, in this approxima- 
tion, the magnetization is linear in A, namely 



A = ^l-JnM 



(11) 



where k = V]£ll. jg the small magnetization limit of 
k(M). The Eq.(|ll|) is a solution if mpfj > IBtt (see 
Eq.(0)). Substituting M from Eq.([Tl|) in Eq.(||), one ar- 
rives at an equation for the gap. This equation can be 
solved in a standard way and the solution is 
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FIG. 2. The zero temperature momentum distribution n, 
for spin up fermions, as a function oi q = and t — cos 6*. 




FIG. 3. The zero temperature momentum distribution n, 
for spin down fermions, as a function of g = ^ and t = cos d. 



The existence of the Fermi surface explains the hnear 
dependence of the specific heat at low temperature: 
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Here N{0) is the density of states on the Fermi surface. 
One can rewrite the 7 = ^ constant in terms of Elliptic 
Integral of the second kind E{a,x) 
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where s < 1 (see Eq.(g|)). Eq.(|l5[) shows that for 
k{M) >> 1 the specific heat constant 7 is small in f- 
superconducting phase, which closely matches the exper- 
iments with ZrZni and URhGe. 

An important experimental fact is that ZrZni and 
URhGe are superconductors at ambient pressure as op- 
posed to the existence of a quantum phase transition 
in UGe2. To comprehend this difference one considers 
the potential @). The quantum phase transition re- 
sults from the existence of a momentum cutoff A, above 
which the potential is repulsive. In turn, the cutoff ex- 
cistence follows from the relation (3 = > 1, which 
is true when the spin-wave approximation expression for 
the spin stiffness constant p = Mpo is used. The spin 
wave approximation correctly describes systems with a 
large magnetization, for example UGe2- But in order to 
study systems with small magnetization, one has to ac- 
count for the magnon-magnon interaction which changes 
the small magnetization asymptotic of p, p = Af^+"poj 
where a > 0. Then for a small M /3 < 1, and the 
potential is attractive for all momenta. Hence for sys- 
tems which, at ambient pressure, are close to quantum 
critical point, as ZrZn^ and URhGe, the magnon self- 
interaction renormalizes the spin fluctuations parameters 
so that the magnons dominate the pair formation and 
quantum phase transition can not be observed. But if 
one applies an external magnetic field, the magnon opens 
a gap proportional to the magnetic field. Increasing the 
magnetic field the paramagnon domination leads to first 
order quantum phase transition. 

The proposed model of ferromagnetic superconductiv- 
ity differs from the models discussed in |^-|l^ in many 
aspects. First, the superconductivity is due to the ex- 
change of magnons, and the model describes in an unified 
way the superconductivity in UGe2, ZrZn^ and URhGe. 
Second, the solution Eq.(pT|) shows that magnetization 
and superconductivity disappear simultaneously. It re- 
sults from the equation of magnetization, which in turn 
is added to ensure that the fermions which form Cooper 
pairs are the same as those responsible for spontaneous 
magnetization. Hence, the fundamental assumption that 
superconductivity and ferromagnetism are caused by the 
same electrons leads to the experimentally observable 



fact that the quantum phase transition is a transition to 
paramagnetic phase without superconductivity. Third, 
the paramagnons have pair-breaking effect. So, the un- 
derstanding the mechanism of paramagnon suppression 
is crucial in the search for the ferromagnetic supercon- 
ductivity with higher critical temperature. For example, 
one can build such a bilayer compound that the spins 
in the two layers are oriented in two non-coUinear direc- 
tions, and the net ferromagnetic moment is nonzero. The 
paramagnon in this phase is totally suppressed and the 
low lying excitations consist of magnons and additional 
spin wave modes with linear dispersion e(fc) ~ fc [ p^ . 
If the new spin-waves are pair breaking, their effect is 
weaker than those of the paramagnons, and hence the 
superconducting critical temperature should be higher. 
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